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Abstract 

We develop a model of ours fitting the quark and lepton masses and mixing angles 
by removing from the model a Higgs field previously introduced to organise a large 
atmospheric mixing angle for neutrino oscillations. Due to the off-diagonal elements 
dominating in the seesaw neutrino mass matrix the large atmospheric mixing angle 
comes essentially by itself. It turns out that we have now only five adjustable Higgs 
field vacuum expectation values needed to fit all the masses and mixings order of 
magnitudewise taking into account the renormalisation group runnings in all sectors. 
The CHOOZ angle comes out close to the experimental bound. 
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1 Introduction 



In previous models |jT|, 3 based on a series of approximately conserved quantum numbers 
taken to be gauged quantum numbers broken by Higgs fields with "small" expectation 
values (VEVs) we managed to fit to typically better than a factor two - on the average 
±60% - all the masses and mixing angles of quarks and leptons known so far. This we 
did by means of six adjustable Higgs VEVs (not counting the Weinberg-Salam Higgs, 
which is fixed from the Fermi constant) assuming that all coupling constants are of order 
of magnitude unity, only Higgs fields having VEVs of orders of magnitude different from 
unity at Planck scale. The adjustable Higgs field VEVs were put up to break a series 
of gauge symmetries - it is enough for the model prediction purpose to consider abelian 
gauge fields, actually two for each family of quarks and leptons - and one of them, 
fits the overall seesaw neutrino scale, while two, p and a;, connect the first and second 
family by breaking the to these famihes associated four ^7(1) gauge groups down to only 
two. Now, however, there were in the previous version of our model three Higgs fields, W , 
T and which in a similar way broke two linear combinations of the U{1) gauge groups 
associated with the families number 2 and 3. But that is seemingly one field too many. It 
is the main purpose of the present article to point out that the model can be improved by 
removing from existence one of those, last mentioned three fields, namely the one called 
X- For making this avoidance of the field, Xi to fit we only need to make an adjustment 
of the quantum numbers of the seesaw scale giving Higgs field (^b-l, which we in the new 
version call 0ss- 

Our newest modification of removing the field, Xi is really only of essential significance 
for the neutrino oscillation parameters because this field only occurred in very insignificant 
matrix elements for the mass matrices of the charged quarks and leptons. By fitting now 
to the neutrinos with fewer parameters in a different way there will of course come a 
slightly different pull for the parameters than before also in the charged sector, though. 

2 The model 

The backbone of the model is the assignment of the in Table 1 described quantum numbers 
under the group 

SMG X U{lf , (1) 

supposed to be conceived of as the gauge group assigning two U{1) gauge groups to each 
family of quarks and leptons. In our scheme with three seesaw neutrinos, a family consists 
of a usual Standard Model family enriched by one right-handed neutrino (to be used as 
seesaw neutrino). For each family the one of these U{1) gauge groups couples just like 
to the weak hypercharge, i.e., with y/2 coupling, but only to the family in question. The 
other one couples to the B — L = "baryon number"- "lepton number" again only for the 
family in question. That is to say we speculate of there being altogether six U{1) gauge 
fields, but now one of them - actually the diagonal subgroup of the three of them - is the 
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weak hypercharge gauge group for the whole system and thus we only have added five 
such groups. 

The model can quite naturally be embellished by the addition to the group of two extra 
SU (2) and two extra SU (3) making the whole Standard Model group including a B — L 
charge, which is gauged, become replicated into one copy for each family of quarks and 
leptons. So each family gets its own gauge fields of all the types in the Standard Model 
and in addition a gauged B — L charge also separate for each family. However, we do 
not need this embellishment for the mass protections that provide the small hierarchies, 
which we fit. For that purpose the abelian quantum numbers in Table 1 are sufficient. 
The embellished model has the gauge group, 

,^^,,3(^MG,xf/(l),_,,0 , (2) 

which we call the family replicated gauge group. 

The quantum numbers of the fermions and of the Higgs fields which we have chosen 
to have in our by now a bit simplified model - because of the deletion of one Higgs field 
- are listed in Table 1. Here the charges, ?/i/2, are the weak hypercharges for the three 
families i = 1,2, 3, while the Baryon number minus lepton number charges B — L for the 
three families separately are denoted {B — L)i, also with i = 1,2, 3. 

The quantum numbers for the embellished model can if one likes be constructed rather 
easily by means of the relation between the non-abelian representation, i.e., of SU{2) and 
SU{3), and the weak hypercharge y/2m the Standard Model. One namely simply requires 
to have the same relation valid for each value of the family denoting index i separately. 
That is to say, that the representation of SU{2)i and SU{3)i for an arbitrary particle, 
quark, lepton, or Higgs, are simply postulated in our embellished version of the model to 
be the same as that quark or lepton Weyl field in the Standard Model that has the value 
y/2 = yi/2 ioY its weak hypercharge (actually, though, we count it mod 1 in the case 0ss)- 

This embellished model has the beauty of having the biggest gauge groups under some 
restrictions such as: only transforming the known particles and the seesaw neutrinos, not 
unifying Standard Model irreducible representations, and having no anomalies 0, |^ . 

2.1 Mass matrices 

With the system of quantum numbers in Table 1 one can easily evaluate, for a given 
mass matrix element, the numbers of Higgs field VEVs of the different types needed to 
perform the transition between the corresponding left- and right-handed Weyl fields. The 
products of Higgs fields needed to achieve the quantum number transitions between left- 
and right-handed fermions are after removal of the Higgs field in the previous version 
of our model, called completely unique. We do no longer have several possibilities 
of achieving the same transition so that a choice of the one giving the dominant matrix 
element would be needed. Rather it is simple linear algebra with the charges to com- 
pute the needed powers for our different fields. The main point of the present type of 
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Table 1: All U{1) quantum charges in the family replicated model. The symbols for the 
fermions shall be considered to mean "proto" -particles. Non-abelian representations are 
given by a rule from the abelian ones (see section 2). 
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model is that the order of magnitudes of the mass matrix elements are determined by 
the number and order of magnitudes of the Higgs field vacuum expectation value fac- 
tors needed. We namely take all the Yukawa couphngs (and also other couplings) in 
our model and the masses of particles that are not mass-protected to be of order unity 
in some fundamental scale (which we usually imagine to be the Planck scale, and actu- 
ally take to be that in the rcnormalisation group runnings of the effective couplings in 
our model). Thereby of course the whole calculation we can do will only give order of 
magnitudewise results! Thus one shall imagine that every expression for a mass matrix 
element in addition to the Higgs vacuum expectation value factors should be provided 
with an unknown factor of order unity (which we in the calculation take to be a ran- 
dom number of order unity), then later logarithmically averaging the results. We have, 
however, not to complicate the expressions too much left out as to be understood these 



3 



factors Aj-J , where q e {up, down, charged lepton, Dirac neutrino, Majorana neutrino}, 
of order unity. The order unity factors physically represent that all the enormously many 
couplings and masses of not mass protected particles are realistically unknown so that 
the only realistic assumption about them is that they are of order of magnitude of unity 
(in fundamental units). In our numerical evaluation of the predictions of our model we 
take the "order unity factors" just in front of each matrix element in our mass matrices 
as dummy integration variables. In fact we integrate the logarithm of, say, a predicted 
mass or mixing angle (as at first a function of our 9 ■ 4 + 6 = 40 order of unity parameters) 
over the real and imaginary parts of these parameters with a weight which is normalised 
so that the corresponding integral of one would be unity. The weight is to be chosen so 
that it emphasises the "order of one factor" being indeed of order unity. This integral 
is then used as our prediction for the logarithm of the mass or mixing in question. This 
integration over the real and imaginary parts of the "order unity factors" is in practice 
performed by a Monte Carlo method. Had we simply put the order one factor equal to 
unity, we would in some simple cases have got roughly right results, however, in many 
cases subdeterminants might become exactly zero by such procedure and we would have 
got severely wrong results - even order of magnitudewise - compared to our physical 
picture of the fundamental couplings being in reality not exactly unity but only of such 
order. 

With the order of unity factor removed our mass matrices look as follows: 
the up-type quarks: 




1 WW 



(3) 



the down-type quarks: 



Mo ~ 



(0ws) 
V2 



\uj^pW\T^Y 



ujp'^WiT^f upW 
I1/2(Tt)4 WT 




the charged leptons: 




(cut)Vl^(Tt)2 

iy(Tt)2 

(iyt)27^4 



(cut)VW^^(T^)' 
WT 




the Dirac neutrinos: 




{u^fp^WT^ 
WW 
{W^)W 



{u^fp^w\T^y 
w\T^y 

wW 





f {p^fT^ (pt)3r6 [p'^fW^T^f 
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The philosophy of the model is that these mass matrices correspond to effective Yukawa 
couplings to be identified with running Yukawa couplings at the fundamental /Planck scale 
for the Higgs field 0ws iii the case of the first three mass matrices and for in the case 
of the right-handed neutrino mass matrix. Therefore these effective Yukawa couplings 
have in principle to be run down by the beta-functions to the scale of observation, see 
section ^ It is also important that we include the "running" of the irrelevant operator of 
dimension 5 giving the neutrino oscillation masses. 

The right-handed neutrino couplings - or mass matrix - is used for producing an 
effective mass matrix for the left-handed neutrinos which we after the five dimensional 
running down mentioned take as the ones observed in neutrino oscillations 0, |^ : 

M,^^M!^M],\M^f. (8) 



2.2 Mcff crude calculation 

Taking into account the various orders of magnitude of our fitting parameters, the ex- 
pectation values in vacuum of Higgs fields, it is not difficult to perform the calculations 
crudely and very often only one or two terms dominate a given quantity. Even the slightly 
more complicated case of the effective matrix for the left-handed neutrinos - the neutrino 
oscillation mass matrix - we compute crudely: 



M{f part 



eff 



A{<Pss) W^T^p^ 



I 



W^p^ 



\ V2W'^T'p^ 



p'WT^ 

W'^p^ 
W^p^ 
^W^T^p^ 



W^T^ WT 

V2W^T^p^ 
V2W^T^p^ 



X 



Lu^p^WT^ 



p^WT^ 



p^W^T^ ' 
WT 



M 



W^T^p^ 



part 
LU'^W'T^ 

V2iu^p^W^T^ 
V2u^p^W^T^ 



V2UJ^P^W^T^ V2cu3^3yf/5y6 

2p^W^T^ 2p^W^T^ 
2p^W^T^ ^p'^W^T^ 



part 



(9) 



Here we used the following rules of calculation aimed at taking into account the random 
numbers of order unity not written explicitly, but meant to be conceived of as standing 
as factors in front of each term in the mass matrices: 



1. When p terms with the same set of "VEV-factors" , i.e., of the types W, T, cj, p 
{4>ws and 0ss) we imagine random phases in the addition (in the complex plane) and 
take the result as y/p times the term that were p times repeated. 

2. We only take the dominant or couple of dominant terms with respect to the VEV- 
factors. 
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3 Renormalisation group equations 



From the Planck scale down to the seesaw scale or rather from where our gauge group is 
broken down to SMG x U{1)b-l we use the one-loop renormalisation group running of 
the Yukawa coupling constant matrices, F^/, Yd, Yg, and (being proportional to the 
mass matrices My, M^,, M^, and M^, respectively), and the gauge couplings [H: 



^dgi 41 3 



^dg2 19 3 

'^"^ = -16^^' ^''^ 

16vr^^ = -7^73^ , (12) 

16vr^^ = I {y,{Y,)' - Yo{Y,)^) Y, + [y^ - (^^? + \gl + 8gl) } Y, , (13) 

le^r^^ = ^ (i^.(r,)^ - r,(r,)t) + { n - (1^? + ^^.^ + sgl)] y, , (i4) 

levr^^ = ^ (l^.(V'.)^ - Y,XY^y) + { n - (^^? + ^^,^)} , (15) 

(K(y;)t - Y,{Y,y) K + { - (^^^^ + ^^7,^)} k , (16) 



IGtt^ 



dt 2 



= ((y;)ty;) + ((K)^^)'^ , (17) 



dt 

Ys = Tr(3FjFy + 3r^ty^ + yjy^ + ytyj ^ (Ig) 

where t = In /i and /i is the renormalisation point. 

However, below the seesaw scale the right-handed neutrino are no more relevant and 
the Dirac neutrino terms in the renormalisation group equations should be removed, and 
also the Dirac neutrino Yukawa couplings themselves are not accessible anymore. That 
means that, from the seesaw scale down to the weak scale, the only leptonic Yukawa 
bet a- functions should be changed as follows: 

16vr^^ = I {YE{Y,y) Y, + [y,- Qgf + ^gj) } Y, . (19) 

Note that the quantity, Yg, must be also changed below the seesaw scale: 

Ys = Tt{3Y^Y, + 3Y^Y, + Y^Y,) . (20) 



Really we stopped the running down according to formula (|I6|) differently for the 
different matrix elements in the Y^, matrix corresponding to the right-handed neutrino 
mass supposed most important for the matrix element in question. 

Starting the running in an analogous way, we further should evolve the effective neu- 
trino mass matrix considered as a five dimensional non-renormalisable term HTH from the 
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different right-handed neutrino masses to the wealc scaleQ (180 GeV) depending on the 
terms: 

levr^^^ = {-"igl + 2A + 2Ys) M,s - ^ {m,s (Y^Y^) + {Y^Y^f Mes) , (21) 

where Yg defined in Eq. (^Of) and in this energy range the Higgs self-couphng constant 
running equation is 

16n'^ = UX' - (^Igl + 9gl) A + ^ (|^^ + ^glgj + gt) + ^YsX - AHs , (22) 

with 

Hs = Tr js {Y^Y.f + 3 {Y^f + {y^Y,)'] . (23) 

The mass of the Standard Model Higgs boson is given by = A {(pws)'^ and, for 
definiteness, we take Mh = 115 GeV at weak scale. 

From 180 GeV down to 1 GeV - experimental scale^ (1 GeV) - we have evaluated the 
beta-functions with only the gauge coupling constans. In order to run the renormalisation 
group equations, we use the following initial values: 

Uil) : giiMz) = 0.462 , giiMp,^^,^) = 0.614 , (24) 
5f/(2) : (?2(Mz) = 0.651 , (72(Mpianck) = 0.504 , (25) 
SU{3) : g^{Mz) = 1.22 , gsiMp,^^,^) = 0.491 . (26) 

Note that we have ignored the influence of the B — L gauge coupling constants; however, 
this effect should not be significant, i.e., Planck scale to the seesaw scale (~ 10^^ GeV) is 
only 10^ order of magnitude difference. Therefore, it should be good enough for our order 
magnitude calculations. 



4 Neutrino oscillations 

The Sudbury Neutrino Observatory (SNO) collaboration has reported p| recently the 
measurement of the neutral current of the active ®B solar neutrino flux and related it to 
measurements of the day and night solar neutrino energy spectra and rates. Moreover, 
they presented improved determinations of the charged current and neutrino-electron 
scattering rate. It turns out that, global analyses of solar neutrino data - combination 
of the SNO results with previous measurements from other experiments |Ty, 0, |1^, 0, |1^, 
15|, |T6| - have confirmed that the Large Mixing Angle MSW (LMA-MSW) solution 



gives the best fit to the data and that the LOW solution is allowed now only at 2.5cr and 



^We take the weak scale as 180 GeV, for simplicity. At this scale we calculated the pole mass of top 
quark, too, using Mt = mt{M) {l + | ^ ■ 

^We set the experimentally observable scale as 1 GeV, thus the charged fermion masses and mixing 
angles are compared to "measurements" at this scale, except the top pole mass. 
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Small Mixing Angle MSW (SMA-MSW) solution at 3.7a, respectively. Not only that we 
"know" the solution of the solar neutrino puzzle but also that these recent results tell us 
that a non-zero CHOOZ mixing angle [|1^] is strongly favoured at the 3.3cr CL. for the 
LMA-MSW solution 



The best fit values of the mass squared difference and mixing angle parameters in the 
two flavour LMA-MSW solution somehow depends on the analysis method, but we take 
the following point as fit values (see in Table 2): 

Am^ = 5.0 X 10"^ eV^ and tan2^0 = O.34. (27) 

The atmospheric neutrino parameters are the following according to the Super-Kamio- 
kande results pO| : 

Aml^^ = 2.5 X 10"^ eV^ and tan^ O^tm = 1-0 . (28) 



5 Results 

The calculation using random numbers and performed numerically was used to fit the 
masses and mixing angles to the phenomenological estimates by minimising what we call 
"goodness of fit" , 

^TT'i ,pred 



In 



(29) 



a kind of for fhe case that we have only order of magnitude accuracy. The result of 
the fitting is presented in Table 2. The results presented there were obtained from the 
following values of the set of Higgs VEVs - where the Higgs field VEVs for the fields p, 
u, T and W causing the breaking to the diagonal subgroup SMG x U{1)b~l are quoted 
with the VEV in "fundamental units", while they are for (p^s (and the not fitted 0ws) 
given in GeV units: 

{W) = 0.157 , (T) = 0.0766 , {to) = 0.244 , (p) = 0.265 , 

(0ss) = 5.25 X 10^^ GeV , (0ws) = 246 GeV . (30) 



The results of the best fit, with the VEVs in Eq. (pO]), are shown in Table 2 and the fit 
has g.o.f. = 3.63. To see a typical error, say average error, compared to the experimental 
values we should divide this value with the number of predictions (17 — 5 = 12) and 



then take the square root of it: y 3.63/12 = 0.55. This means that the 12 degrees of 
freedom have each of them a logarithmic deviation of about 55%, i.e., we have fitted all 
quantities with a typical error of a factor exp (^3.63/12) ^ 1.73 up or down. This agrees 



with theoretically predicted deviations [^]. However, our worst fitting value is now the 
electron mass which we predict/fit a factor 3 too heavy. 
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Table 2: Best fit to conventional experimental data. All masses are running masses at 
1 GeV except the top quark mass which is the pole mass. Note that we use the square roots 
of the neutrino data in this Table, as the fitted neutrino mass and mixing parameters, in 
our goodness of fit (g.o.f.) definition, Eq. (PU]). 
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We should here emphasise this point: even though we reduced the number of free 
parameters by one - now only five VEVs - we are able to fit all fermion masses and 
mixing angles with an average factor 1.73 of deviation as in previous work. 

Experimental results on the values of neutrino mixing angles are usually presented in 
terms of the function sin^ 29 rather than tan^ 9 (which, contrary to sin^ 29, does not have 
a maximum at 9 = n/A and thus still varies in this region). Transforming from tan^ ^ 
variables to sin^ 29 variables, our predictions for the neutrino mixing angles become: 



sin2 2^0 = 0.66 , (31) 

sin2 2^atm = 0.96 , (32) 

sin2 2^ehooz = 0.11 . (33) 

We also give here our predicted hierarchical left-handed neutrino masses (mj) and the 
right-handed neutrino masses (Mj) with mass eigenstate indices [i = 1,2,3): 

nil = 1.4 X 10-3 eV , Mi = 1.0 x 10^ GeV , (34) 

ni2 = 9.6 X 10-3 eV , Ms = 6.1 x 10^ GeV , (35) 
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ms = 4.2 X 10-2 eV , M3 = 7.8 x 10^ GeV . (36) 

Note that our fit of the CHOOZ angel (Eq. ^3]) is lying on the borderline of the 
experimental results which analysed by two flavour method. However, our fit satisfy even 
2(T C.L. limit (tan^ 6'chooz = 3.3 x 10"^) based on three flavour analysis ||2^ . 



For the CP- violation parameter, the Jarlskog triangle area [^], Jcp, we expect our 
prediction to have a larger uncertainty than for the quantities, which like the masses are 
essentially mass matrix elements, because it corresponds to a ratio or product of six such 
quantities and thus our prediction, 

JcP,p.ed = 3.6 X 10-^ (37) 
compared to the "experimental" value 

^cP.o.p = (2 - 3.5) X 10-^ , (38) 

should be considered a success of only l.So" deviation. 

For the for the neutrinoless double beta-decay relevant specially weighted neutrino 
mass average 

3 



=1 



(39) 



we obtain the prediction from our fit 

|(m)| = 3.5 X 10-^ eV . (40) 

Our proton decay predictions may vary a bit with the philosophy, however, we in any 
case predict so long lived protons that it is pretty hopeless to see any proton decay. At 
least we get a proton life time of the order: 

^(p^^Oe+) ~ 10^^ years . (41) 



The brenching ratio of — e + 7 is performed by |23] 



where My/ is the mass of the bosons. We insert our predictions from Table 2 and 
Eq. ( ^11) in this formula, then we get 

Br(/i ^ e + 7) ~ 10-^^ . (43) 
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6 Simple relations 



In the approximation of only taking dominant terms and imagining all the quantities 
extrapolated to the Planck scale - i.e. ignoring renormalisation group running and sub- 
dominant terms - we can get the rather simple relations: 

1. Family degeneracy for charged quarks and leptons: 

rrib ~ rrir , ^ , rriu ^ ^ rrie , (44) 
(However, we avoid the rric and rrit being degenerate with their families). 

2. Factorisation of mixing angles for quarks: 

VcbVus^V^b , (45) 

3. Neutrino relations: 

^O^^atm ^ ^^coo. , (^^) ^ , ^ , (46) 

Here mi is the mass of the lightest left-handed neutrino. 

4. Relations for charged quarks and leptons: 

2 

ml ^ mt nis nic , rrit ^ "Weak scale" , Vcb ~ — ■ (47) 

mcrrib 

5. Relations between neutrinos and charged quarks and leptons: 

^atm ^ V2— ^ V2—^^ , e^^-^ — ^ . (48) 



Note that according to Eq. (^61) the neutrino data fit by just by two parameters, ^0 
and 9a.tm, and one overall scale, (0ss)- Moreover, according to Eq. ( |48| ) the mixing angles 
are given in term of the quark quantities. 



7 Theoretical lesson for the mass matrix 

We should emphasise that we have obtained the rather good fit for all the neutrino 
quantities - 9q, 6'atm, 6'chooz, AttIq and Am^^^ - and the charged masses and mixing 
angles using only 5 free parameters, namely, five VEVs of Higgs fields. Our effective 
left-handed neutrino mass matrix, Eq. (H), is of the form 

^ (pi 0102 0103 \ 

0102 0i 0203 , (49) 

V 0103 0203 03 / 
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with different order unity factors, wliicli we call factorised mass matrix. It is as a con- 
sequence of tliis type of form we get a relation (|46|) . In fact, the ratio of mass squared 
difference, AmQ/Am^^^^, is characterised by the fourth power of the atmospheric mixing 
angle at Planck scale. Due to the large atmospheric neutrino mixing angles, it is clear 
that this fact - mass squared difference being order unity - is not true experimentally 
if the effects of renormalisation equations are not taken into account on the Dirac- and 
Majorana-sectors from Planck scale to seesaw scale, and also non-renomalisable five di- 
mensional evolution equation from the seesaw scale to the weak scale. However, with 
the running included the "factorised" mass matrix is a good candidate for fitting the 
LMA-MSW solution due to only moderate hierarchy structure of this solution. 



8 Baryogenesis not so successful 



Our prediction of baryogenesis in the Fukugita-Yanagida scheme p5| is unfortunately not 
so successful: The asymmetry from the decay of the lightest of the seesaw neutrinos is 
way too small to produce the B — L to give the baryon number asymmetry needed for 
the light element production in big bang at minute scale. However, the asymmetry from 
the decay of the two heavy and approximately degenerate seesaw neutrinos is much larger 
and so our hope would be that they could produce the asymmetry needed; however, it 
turns out that the wash-out by the processes involving the lightest right-handed neutrino 
becomes too large and we get too little baryon number at the end because of the rather 
high value of the mi [p6| . 

It must though be said that the m's are products or ratios of three mass matrix ele- 
ments so that the uncertainty is expected to be exp(±A/3 ■ 64%) = exp(±lll%), meaning 
a factor 3 up or down. Our problem with the rather high value of m's is partly due to 
the renormalisation group corrections without which lower m's by about a factor 3 would 
result. We also hope that taking into account that different flavours will be washed out 
with different rates will lower the wash-out rate at the end. Somewhat optimistically 



this could mean that we could use an effective mi a factor 3 smaller. 

Taking into account also the uncertainty in other quantities we might hope for that 
the baryon number predicted only disagrees by about three standard deviations. 

The many questions of discussion involved in calculating the baryon number we should 
say that we got pretty tight to a viable result, however, only based on the idea that it is 
the heavy seesaw neutrinos that produced the baryon asymmetry. 



9 Conclusion 

We have developed our previous model by deleting from the list of assumed Higgs fields 
with vacuum expectation values breaking the assumed gauge group SMG x f/(l)^ or 
iJi,2,3{i3MGi X a Higgs field, X; which were already suspected to be unwanted 
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by having the same quantum numbers as the combination W^(T^)^. We found that we 
can indeed fit to order of magnitude accuracy the nine charged quark and lepton masses, 
the two neutrino mass square differences, the five measured mixing angles in addition 
to matching the CP-violation with only five Higgs field vacuum expectation values being 
adjusted to order of magnitude accuracy. In addition this fit manages to pass the test 
of several experimental bounds: no neutrinoless double beta-decay, no proton decay, and 
no n —>■ e + 'J within present experimental accuracy, and most interestingly we predict a 
sufficiently small CHOOZ angle, U^e] however, there is the prediction that we are actually 
close to the limit so that our model might be realistically falsified by decreasing the 
CHOOZ angle limit by an order of magnitude. 

It is remarkable that we fit such a large number, 17, of genuinely measured quantities 
with only our 5 parameters order of magnitudewise. Our precise gauge group proposals 
are not uniquely called for in as far as we could at least include more or less embellishment 
with non-abelian groups for the different families or replace some of the abelian groups by 
non-abelian, but the idea that different families have different "chiral" quantum numbers 
may be hard to avoid without throwing our good agreement out as being accidental. 
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